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Connections between charge balance functions, charge fluctuations and correlations are presented. 
It is shown that charge fluctuations can be directly expressed in terms of a balance functions under 
certain assumptions. The distortion of charge balance functions due to experimental acceptance is 
discussed and the effects of identical boson interference is illustrated with a simple model. 



I. INTRODUCTION 



II. RELATING BALANCE FUNCTIONS, 
FLUCTUATIONS AND CORRELATIONS 



Charge balance functions and charge fluctuations 
H ||, ||, U have been proposed as a means for gaining 
insight into the dynamics of hadronization in relativis- 
tic heavy ion collisions. Both observables are sensitive to 
the separation, in momentum space, of balancing charges. 
Such a pair is composed of a positive and negative par- 
ticle whose charge derives from the same point in space- 
time. As a quark-gluon plasma scenario entails a large 
production of new charges late in the reaction, a tight cor- 
relation between the balancing charge/anti-charge pairs 
would provide evidence of the creation of a novel state of 
matter. 



Both balance functions and charge fluctuations can be 
expressed in terms of one-particle and two-particle ob- 
servables. In the following section we present expres- 
sions for both balance functions and charge fluctuations 
in terms of spectra and correlations, and show how the 
charge fluctuations can be simply expressed in terms of 
balance functions for neutral systems. 

Unfortunately, balance functions and charge fluctua- 
tions can both be rather sensitive to detector acceptance. 
In Sec. 



[II we present a variant of balance functions 
which reduces acceptance effects for detectors with sharp 
cutoffs in rapidity. Identical pion correlations also af- 
fects both observables in a non-trivial manner. In a sim- 
ple model utilizing parameters consistent with observed 
spectra and correlations from RHIC, we illustrate the 
distortion of the balance functions due to Bose-Einstein 
correlations in Sec. IV. The insights gained from these 
studies are summarized in Sec. |v|. 



* Electronic address: 
t Electronic address: 



jeon@hcp.phvsics.rr] 


cgill.ca 


3ratt@nscl.msu.edu 





As mentioned in the introduction, both balance- 
function and charge-fluctuation observables are gener- 
ated from one-body and two-body observables which ne- 
cessitates that they may be expressed in terms of spectra 
and two-particle correlation functions. In order to ex- 
press the balance functions in terms of the elementary 
correlation functions, first define 
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where A 1,2 are phase space criteria such as rapidity in- 
tervals. 

In terms of these quantities, the balance function is 
expressed as 

S(A 2 |Ai) = i{Z>(- A a |+, AO -£>(+, A a |+,Ai)(5) 
+ £>(+, A a |-,Ai) -£>(-, A 2 |- A x )}. 



where 
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which can be considered as a conditional probability. 

Thus correlation functions and spectra are sufficient 
to determine balance functions, although the required 
integration can be somewhat convoluted, depending on 
the binning, Ai and A 2 . The criteria Ai is based solely 
on the momenta of the first particle, while the criteria A 2 
might be any function of the momenta of both particles, 
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e.g., it might be determined by the relative rapidity of 
the two particles. 

To establish the correspondence between charge fluc- 
tuations and balance functions, consider a balance func- 
tion binned as a function of the rapidity difference where 
both particles are required to reside within a fixed ra- 
pidity window of size Y. For this case Ai constrains the 
first particle to be within the rapidity window, and A2 
constrains the second particle to have a relative rapid- 
ity \yb — y a \ = Ay while also existing inside the rapidity 
window. This binning was applied in preliminary results 
from STAR reported in ||. Referring to this balance 
function as B(Ay\Y), one can find the charge fluctua- 
tion within the rapidity window < y < Y by integrating 
B(Ay\Y) in the interval < Ay < Y. In this case, 



B(Y\Y) = 



dAy B(Ay\Y) 
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where (■ ■ -)a denotes averages in the phase space region 
A. Writing N± = (N±) A + 5N±, it is not hard to show 
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where Q = N + - and N ch = N + + N_. For electric 
charge, the size of the correction is usually less than 5 % 
in relativistic heavy ion collisions where the number of 
produced charges is much greater than the net charge. 
However, for baryon number the additional term is not 
negligible even at RHIC. 

In a boost-invariant system (independent of rapidity) 
the balance function B(Ay\Y) can be related to the bal- 
ance function for an infinite interval. 



B(Ay\Y) = B(Ay\Y = oo)(l - Ay/Y). 



(9) 



The factor (1 — Ay/Y) accounts for the probability that 
a particle's partner will fall within the rapidity window 
given that they are separated by Ay. Also, assuming 
boost invariance allows one to express the balance func- 
tions simply in terms of correlation functions as described 
in Eq. (§. 
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From the above discussion it is clear that the charge 
fluctuation is the global measure of the charge correla- 
tion and the balance function is a differential measure 
of the charge correlation and therefore carries more in- 
formation. The advantage of charge fluctuations is that 



they carry a clear physical meaning in terms of a grand 
canonical ensemble [|[ ||, and can therefore be easily 
connected to more ideal theoretical models, e.g. Lattice 
QCD calculations. However, since there are no exter- 
nal sources of charge in heavy ion collisions to warrant 
a grand canonical treatment, both observables are effec- 
tively driven by the dynamics of how balancing charges 
are formed and separate. 

We emphasize here that charge fluctuations were not 
intended toprovide a derivative measure. As can be seen 
from Eq. (|3|) the charge fluctuation summarizes the bal- 
ance functions in one number. It gives somewhat differ- 
ent information than the width of the balance function 
since it is also affected by the height. We do not rec- 
ommend analyzing charge fluctuations as a function of 
the size of the rapidity window. If the different sized 
windows included the same pairs, the values would no 
longer be statistically independent when plotted against 
the window size. If the windows are used only once, the 
information from pairs which occupy adjacent windows 
is thrown away. 

A similar set of issues surfaced in making the connec- 
tion between fluctuations and correlations in the study of 
multiplicity distributions analyzed as a function of rapid- 
ity 0, g|. A more general connection between fluctuation 
and inc usive observables can be found in qJ. However, it 
should be noted that factorial moments and scaled facto- 
rial moments, which are measures of fluctuation 0, [ll), 
offer the opportunity to study n-body correlations for 
n > 2 in a manner which, unlike correlations, can be 
easily collapsed into a single variable. 



III. MINIMIZING ACCEPTANCE EFFECTS IN 
BALANCE FUNCTIONS 

Balance functions analyzed by the STAR collabora- 
tion H were constructed according to the prescription 
that pi would refer to any pion that is measured within 
a specified rapidity window while -pi referred to the rela- 
tive rapidity, again with the requirement that the second 
particle was within the rapidity window. In that case, 
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Here (Ay) counts pairs with opposite charge that 

satisfy the criteria that their relative rapidity equals Ay, 
whereas N + is the number of positive particles in the 
same interval. Here the angular bracket represents av- 
eraging over the events and Y is the size of the detec- 
tor rapidity window. From this example, one can read- 
ily understand how balance functions identify balancing 
charges. For any positive charge, there exists only one 
negative particle whose negative charge derived from the 
point at which the positive charge was created. By sub- 
tracting from the numerator the same object created with 
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positive-positive pairs, one is effectively subtracting the 
uncorrelated negatives from the distribution and identi- 
fying the balancing charge on a statistical basis. 

From the construction of Eq. (pd|), one can understand 
the sensitivity of B(Ay\Y) to the acceptance size Y by 
considering a detector which covers a finite range in ra- 
pidity. 
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with y max — 2/min — Y . For this example, the balance 
function must go to zero as Ay approaches Y. This oc- 
curs because the particle satisfying the condition p\ must 
lie at the extreme boundary of the acceptance in order 
for the second particle to have a relative rapidity Ay ~ Y 
while remaining in the acceptance. The balance function 
is thus forced to zero at the limits of the acceptance for 
trivial reasons. 

Of course, the balance function corresponding to a per- 
fect detector, B(Ay\oo), is independent of any particu- 
lar detector size Y. As described in Eq. (||), one can 
easily correct for the detector acceptance in the boost- 
invariant case by dividing the balance function by a fac- 
tor (1 - Ay/Y). 

These balance functions would not have more informa- 
tion than those created without the correction factor, but 
the information would more directly address the physics 
of charge separation rather than reflecting the experimen- 
tal acceptance. We note that the statistical uncertainties 
of the corrected balance function will however be quite 
large as Ay ~ Y. 

More generally, one can correct the balance func- 
tions for the acceptance by dividing the numerators 
in Eq. (||), N(Q 2 ,p 2 \Qi,pi), by acceptance factors, 
A(Qi,p2\Qi,pi). 
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The acceptance factor represents the probability that, 
given a particle i satisfied the criteria p\ , a second parti- 
cle that satisfied p 2 would be detected. Since the criteria 
(Q 2 ,p 2 ) may depend on the individual particle that satis- 
fied (Qi,Px), it may be simpler to calculate A in terms of 
o-i(Q 2 TP 2 ) which represents the acceptance into (Q 2l p 2 ) 
given the particular particle i. 

EieQ 1)Pl a l {Q 2 ,P2) 



A(Q 2 ,p 2 \Qi,pi) 



(14) 
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The acceptance probability di(Q 2 ,p 2 ) would be between 
zero and unity. We note that the acceptance is effec- 
tively accounted for by performing a substitution for the 
denominators in Eq. (|5|). 

AT(Qi,i\)-> Yl <Q^) (15) 

i£Ql,Pl 



For the boost-invariant case above where p 2 referred 
to the relative rapidity, and where the acceptance is rep- 
resented by simple step functions in rapidity, the proba- 
bilities would become 



ai(p 2 ) 




IX 

- yi < Ay and j/j - y min < Ay 

- yi < Ay and y l - y min > Ay 

IX 

- yi > Ay and y l - y min < Ay 
IX - Hi > Ay and y l - y min > Ay 

(16) 

Given that the bins p 2 would be of finite extent, the val- 
ues might differ from 1/2, or unity if the bin straddled 
the acceptance. For a boost invariant system, averaging 
over yi results in the simple correction factor (1 — Ay/Y) 
mentioned previously. 

In general, if the acceptance depends on where in the 
{Q 2 ,P 2 ) bin the second particle lies, one can not cal- 
culate the acceptance correction exactly without knowl- 
edge of the charge correlation which is unavailable except 
by measuring the balance function in sufficiently small 
p 2 bins such that the acceptance is effectively uniform 
within the small bins. This may not be feasible due 
to statistics. It is our recommendation that such fac- 
tors cii(Q 2 ,p 2 ) should be kept simple. One can always 
correct theoretical results for the detector response by 
applying whatever factor is applied to the experimental 
analysis. Although comparisons with models could have 
been made without any corrections, acceptance-corrected 
balance functions can allow for a more physical interpre- 
tation while not compromising the integrity of the anal- 
ysis. 



IV. BOSE- EINSTEIN CORRELATIONS AND 
BALANCE FUNCTIONS 

Although Bose-Einstein correlations only affect iden- 
tical particles at small relative momentum, they mani- 
fest themselves in balance functions despite the fact that 
the binning in balance functions typically covers a large 
volume in momentum space. In a related topic, Bose- 
Einstein correlations (also known as the Hanbury-Brown 
Twiss effect, HBT) have been observed in rapidity cor- 
relations where all charged particles, both positive and 
negative, were used in the analysis j?], |l2|, |l3| . The man- 
ifestations of HBT in balance function derives from the 
fact that it induces a correlation between a given charge 
and all other charges, not just those that were created to 
balance the given charge. 

In balance functions the HBT effect should enhance the 
probability that same-charge particles have small relative 
momentum, thus providing a dip in the balance function 
at small relative rapidity. In order to model this effect, 
we consider pairs of pions with momenta p a and pi, and 
opposite charge that are created according to a boost- 
invariant thermal distribution with a temperature of 190 
MeV, thus roughly reproducing the pion spectra mea- 
sured in Au + Au collisions at RHIC. In addition to the 
usual contribution to the balance function between p a 
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and pf,, a second component derives from the interaction 
with other pions from other pairs which in this case have 
momenta p c and pd- The thermal distribution describing 
the first two particles was centered at zero rapidity, while 
the thermal distribution responsible for emission of the 
second pair was randomly chosen within ±4 units. 

The particles p a and p c were assumed to have the same 
sign, as were the particles with momenta pb and pd- A 
contribution to the balance function was constructed us- 
ing these particles, but with a weight, 

w = C ++ {p a ,p c )C — (p b ,pd)C + -(p a ,pd)C^+{p b ,p c ). 

This accounts for the weight due to two-particle interac- 
tions. The correlation functions were simple functions of 
Qinv(Pa,Pb) = y/(Pa - Pb) 2 , which were generated by cal- 
culating correlation functions for a spherically symmetric 
Gaussian source of radius, i?i nv = 7 fm, again crudely in 
line with measurements at RHIC[Q. The weights were 
calculated by averaging the squared relative wave func- 
tion for two particles, including the Coulomb interaction 
between the pions. The weight was multiplied by the 
number of such pairs which came from assuming that 
there were 200 pion pairs per unit rapidity. Only a frac- 
tion, A = 0.7, of the pairs were assumed to interact due 
to the fact that some pions would be created in long- 
lived decays. The acceptance of the STAR detector, and 
the fact that only a fraction of the pions would truly be 
balanced by other pions (rather than by charged kaons 
or other particles) was roughly accounted for by accept- 
ing only 60% of the particles with transverse momenta 
between 100 MeV/c and 700 MeV/c. 

The resulting balance functions are displayed in Fig. [I]. 
When the interactions between particles is neglected, the 
resulting balance function falls monotonically, and has a 
width consistent with the temperature. The inclusion of 
the HBT effect results in a large dip near Ay = 0, and 
an enhancement at somewhat larger Ay. The dip derives 
from the enhancement of same-sign pairs which results 
in a negative contribution to the balance function. Since 
the weight is assigned to the emission of the cd pair, the 
positive HBT weight contributes to opposite sign pairs 
with equal strength, but is spread out over a wider range 
of Ay. Hence, the balance function is slightly enhanced 
for Ay ~ 1/2 from HBT effects. 

Also shown in Fig. [l] are calculations where the 
Coulomb interaction is included. Since Coulomb inter- 
actions result in attractions for opposite-sign pairs, and 
repulsions for same sign pairs, the dip due to HBT is 
mitigated. 



Although the shape of the balance function is visibly 
altered by the inclusion of two-particle interactions, the 
mean width changed by only a few percent. The strength 
of the distortion was proportional to the multiplicity, but 
the effect is not necessarily weaker for peripheral events. 
This follows because the HBT correction contributes with 
a strength proportional to R~ 3 . Since the product of 
dn/dy and R~ 3 stays roughly constant over a wide range 
of centralities in heavy ion collisions, the distortion due to 
interactions should not appreciably affect the centrality 
dependence of the balance function's width. 



V. SUMMARY 

This paper covered several technical issues related to 
balance functions. The conclusion of Sec. II is that 
charge fluctuations can be related to balance functions in 
a straightforward manner unless the average net charge 
is large. In fact, the charge fluctuation can be thought of 
as a measure of the integrated balance function B(Ay\Y) 
from zero to Y. In that sense, it represents a one- 
component measure of the balance function, just like the 
mean width. 



Section III provided an illustration of how balance 



functions can be created in such a way as to minimize 
sensitivity to experimental acceptance. Although the ex- 
ample addressed only problems with finite acceptance in 
rapidity where the balance function was binned accord- 
ing to relative rapidity, the principles could be applied to 
balance functions in any variables. 

The last section considered the inclusion of two- 
particle interactions into Balance functions. The effects 
were shown to be quite visible at small relative rapid- 
ity. However, the width of the balance function was 
not significantly affected by the two-particle interactions. 
This is encouraging, as it justifies interpreting balance 
functions as objects that statistically identify balancing 
partners, while subtracting out contributions from other 
pairs. 



Acknowledgments 

This work was supported by the National Science 
Foundation, Grant No. PHY-00-70818, by the Natural 
Sciences and Engineering Council of Canada and by le 
Fonds pour la Formation de Chercheurs et PAide a la 
Recherche du Quebec. 



[1] S. Bass, P. Danielewicz, and S. Pratt, Phys. Rev. Lett. 

85, 2689 (2000). 
[2] S. Jeon and V. Koch, Phys. Rev. Lett. 85, 2076 (2000). 
[3] S. Jeon, V. Koch, K. Redlich, and X. N. Wang, LANL 

report (2001), nucl-th/0105035. 



[4] C. M. Ko et al., Phys. Rev. Lett. 86, 5438 (2001), nucl- 
th/0010004. 

[5] M. Asakawa, U. Heinz, and B. Muller, Phys. Rev. Lett. 

85, 2072 (2000), hep-ph/0003169. 
[6] M. B. Tonjes, Proceedings, 17th Winter Work- 



■5 




Ay 



shop on Nuclear Dynamics, Park City Utah (2001), 

http://theo08.nscl.msu.edu/winterworkshop 

/procecdings/ParkCity _tonjes.htm. 
[7] K. Wieand, S. Pratt, and A. Balantekin, Phys. Lett. 

B274, 7 (1992). 
[8] K. Haglin and D. Seibert, Phys. Lett. B273, 211 (1991). 
[9] A. Bialas and V. Koch, Phys. Lett. B456, 1 (1999), nucl- 

th/9902063. 

[10] A. Bialas and R. Peschanski, Nucl. Phys. B308, 857 
(1988). 

[11] D. Seibert, Phys. Rev. D41, 3381 (1990). 
[12] R. Wilkes and the EMUOl Collab., Nucl. Phys. A544, 
153c (1992). 

[13] Y. Akibaet al., Phys. Rev. C56, 1544 (1997). 
[14] C. Adler et al. (STAR), to appear in Phys. Rev. Lett. 
(2001), nucl-ex/0107008. 



FIG. 1: The balance function from the simple thermal 
Bjorken model (line) has been parameterized and filtered 
to roughly provide rough consistency with measurements of 
STAR. The inclusion of HBT effects (triangles) gives a dip 
at small Ay, while the extra addition of Coulomb (circles) 
modifies the dip. 



